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Abstract 

We apply Painleve test to the most general variable coefficient nonlinear 
Schrodinger (VCNLS) equations as an attempt to identify integrable classes 
and compare our results versus those obtained by the use of other tools like 
group-theoretical approach and the Lax pairs technique of the soliton theory. 
We present explicit transformation formulae that can be used to generate 
new analytic solutions of VCNLS equations from those of the integrable NLS 
equation. 

1 Introduction 

Variable coefficient extensions of nonlinear evolution type equations tend to arise in 
cases when less idealized conditions such as inhomogeneities and variable topographies 
are assumed in their derivation. For example, variable coefficient Korteweg-de Vries 
and Kadomtsev-Petviashvili equations describe the propagation of waves in a fluid 
under the more realistic assumptions including non-uniformness of the depth and 
width, the compressibility of the fluid, the presence of vorticity and others. While 
these conditions lead to variable coefficient equations, all or some of the integrability 
properties of their standard counterparts, namely when the coefficients are set 
equal to constants, are in general destroyed. However, when the coefficients are 
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appropriately related or have some specific form, the generalized equation can still be 
integrable as evidenced by the presence of Painleve property, Lax pairs, symmetries 
and other attributes of integrability [H E]. This fact suggests us to look for a point 
transformation mapping VCNLS equation to the standard NLS one. It is thus our 
intention to examine the integrability conditions of a class of variable coefficient 
nonlinear Schrodinger (VCNLS) equation which can be thought of as a generalization 
of the standard integrable nonlinear Schrodinger (NLS) equation. More precisely, we 
ask the question: Under which conditions can the equation be transformed to its 
constant coefficient version which is notoriously integrable? 
In this paper we consider the following class of equations 

ii!t + fix,t)ip^^ + g{x,t)\i)\'^ip + h{x,t)ip = 0, 
f = fi+if2, 9 = 9i + W2, h = hi + ih2, (1.1) 
f„g„h,eR, J = 1,2, /i^O, gi^O 

and study its integrability and in particular reducibility to its standard form using 
the tools of singularity analysis and symmetry which are well-established approaches 
towards integrability. This class of equations model a large variety of physical 
phenomena and has been widely analyzed through various methods. The complex 
function ijj{x,t) has different physical meanings in different physical contexts. An 
equation with an additional quintic nonlinearity when the coefficient functions mostly 
depend on time models the propagation of pulses in optic fibers and was studied in [3] , 
where trigonometric type solutions were obtained through some transformations. [1] 
studies the equation with space and time coordinates switched, modeling propagation 
of pulses in optical fibers with distributed dispersion and nonlinearity and finds 
soliton-type solutions with a Darboux transformation, whereas elliptic-type solutions 
through various transformations [5], soliton-type solutions via Hirota method [6] and 
in terms of a double- Wronskian determinant [7J, and self-similar solutions also exist 
in the literature [8]. In these works, the coefficients are considered as functions of a 
single variable. For another subclass of equations when /, g are real functions of t, 
reductions to a nonlinear ODE by an appropriate transformation is performed and 
solutions of hyperbohc, trigonometric and elliptic-type are obtained [9]. Besides, 
elliptic-type solutions by a direct-symmetry method are found in [10] for the same 
case, ^llj considers the case with constant /, g and a real potential periodic in time, 
and presents soliton-type solutions. For / = 1 and g, h real, [12] uses similarity 
transformations to transform the equation into the stationary nonlinear Schrodinger 
equation and hence obtains soliton-type solutions of various characters. 

For /, g real function of t and a real potential which is quadratic in x with 
time- dependent coefficients, [13] gives an integrability condition based on a Lax pair. 
[11] tries to transform this equation directly to the NLS equation, whereas [TB] states 
that this equation and its extended version are equivalent to the NLS equation. 

There is a vast amount of literature devoted to analyzing Painleve property 
(P-property) of equations belonging to the class (11. ip as subcases (for example, see 
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[T6| [T7]). Results of [18] and [19] belong to a wider class and the compatibility 
condition they find coincides with the exact integrability condition given by |13j . 
Finally, we would like to mention Ref . [20] , in which allowed transformations of ( ll.ip 
are found and used to do a complete symmetry group analysis of the equation. 
In addition, we apply the Painleve test to the more general equations 

i^t + f^xx + 9\i^\^^ + h^ + kij^ = (1.2) 

and 

i^t + f^xx + 9\^\'^^ + h^ + ktp^ + il^xxx + im{\ip\'^ip)-, + in^{\^\'^)^ = 0, (1.3) 

where the coefficient functions are considered to be arbitrary in the form / = fi{x,t) + 
i f2{x, t), g = gi{x, t) + i g2{x, t), h = hi{x, t) + i h2{x, t), k = ki{x, t) + i k2{x, t), 
I = l{t), m = m{t), n = n{t). P-test results for (11. 2p are obtained in ^21j when 
/ is real. In Ref. |22j integrability of (11.31) is investigated for the coefficients 
having only time dependence. Here we allow also spatial dependence and obtain 
the Painleve-integrability conditions. We refer the interested reader to these works 
for the physical motivation of these generalized equations. 

The paper is organized as follows. Section 2 performs the Painleve analysis of 
equations (ll.ll) - (ll.3p . Section 3 is devoted to the transformation of the equation into 
the NLS equation using allowed transformations. We end up by applying our results 
to a generalized Gross-Pitaevskii equation in a harmonic trap in Section 4. 



2 Singularity Analysis of the Equations 



2.1 Painleve test for eq. ( 11.11 ) 



For convenience we write (ll.ip together with its complex conjugate as the system 

iut + f{x, t)u^^ + g{x, t)u'^v + h{x, t)u = 0, 
—ivt + p{x, t)vxx + Q'(a;, t)uv'^ + r(a;, t)v = 0. 



Here u was employed instead of ip and v denotes its complex conjugate, but in 
this context they are viewed as independent functions, whereas p, q, r are complex 
conjugates of /, g and h, respectively. This system will be subjected to the Painleve 
test for partial differential equations. We shall look for solutions of the system (12. ip in 
the form of a Laurent expansion (known as Painleve expansion) in the neighborhood 
of a non-characteristic, movable singularity surface (actually a curve in this case) 
defined by 

$(x,t) = 0. (2.2) 

Thus, we expand 

oo oo 

u{x,t) = Mj(x,t)$"+^(x,t), v{x,t) = Vj{x,t)<l>'^+\x,t), (2.3) 

3=0 3=0 
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where uo,vo 7^ and Uj,Vj,^{x,t) are analytic functions, a and /3 are negative 
integers to be determined from the leading order analysis. 

The partial differential equation (PDE) is said to pass the Painleve test if the 
substitution of the above expansion into the equation leads to the correct number 
of arbitrary functions as required by the Cauchy-Kovalewski theorem given by the 
expansion coefficients, where ^{x,t) should be one of the arbitrary functions. The 
coefficients in the expansion where the arbitrary functions figure are known as the 
resonances. 

If a PDE passes this test, then it has a good chance of having the Painleve 
property: solutions (12. 3p are single valued about the singularity surface, and depends 
on sufficient number of arbitrary functions which are needed to satisfy arbitrary 
Cauchy data imposed for some x = Xq. We note that passing the Painleve test 
is necessary, but not sufficient for having the Painleve property. In that case, for 
integrable PDEs it is usually possible to construct auto-Backlund transformations 
which relate equations to themselves via differential substitutions and also Lax pairs, 
which then ensures the sufficient condition of integrability. Application of the Painleve 
expansion to nonintegrable PDEs can allow particular explicit solutions to be obtained 
by truncating the expansion which then imposes constraints on the arbitrary functions 
and the function This usually requires compatibility of an overdetermined PDE 
system. Let us mention that the method of truncated expansion has been successfully 
applied to many nonintegrable PDEs in constructing exact solutions. 

For the determination of leading orders a and f3, we substitute u ~ mo$° and 
V ^ Vq^^ in (12. ip and see that by balancing the terms of smallest order 

a + /3 = -2 (2.4) 



must hold, which only allows the negative integers a = —1 and /3 = — 1. In addition, 
we obtain the relations 

Uoi9UoVo + 2f^l) = (2.5a) 
voiquovo + 2p<^l) = (2.5b) 



and it follows that 



^2/ _ n^2P 



uovo = -2^i'- = -2^i-. (2.6) 
9 Q 

From this it is seen that ^ must be real, which requires that fig2 = f29i- We can say 
that, if one of / and g is real or pure imaginary, then so is the other. 

After determination of the leading orders, we substitute (12. 3p into (12. ip . For j = 0, 
equating the coefficient of to zero we exactly obtain the relations (12.50 . For j > 1, 
equating coefficient of $~^+-' to zero yields the system 
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where 

Fj = - iuj-2,t - ^(-2 + i)(j)tUj-i - fuj^2,xx - 2(-2 + j)f^.j,Uj-i^x 

- (-2 + - hUj^2 - gVo ^ UlUj-i - g'^^ UiUk-lVj_k, 

1=1 k=i 1=0 ^2.8) 

Gj = lVj^2,t + «(-2 + j)(t)tVj-l - pVj-2,xx - 2(-2 + j)p^xVj-l,x 

j-l j-l k 

- (-2 + j)p^xxVj^i - rvj-2 - quo ^ t^z^^j-i - viVk-iUj-k- 

1=1 k=l 1=0 

Here we have used the convention that uj = for j < and sums which have upper 
indices that are less than the lower indices will be taken as zero. Notice that Fj and 
Gj can be obtained from each other by the interchange i -ir^ —i, u -H- v, f -ir^ p, g -(r^ q, 
h r. Uj and Vj are determined by the system fl2.7l) unless 

detQ(j) = \f\^^lU + - 3)(j - 4) = 0. (2.9) 

Hence we find the resonance levels to be j = —1,0,3,4. The resonance j = — 1 
corresponds to the arbitrariness of the singularity manifold $ and the resonance 
j = points out that there is one arbitrary function among uq and Vq, as it is verified 
by (S. 

For j = 1 we solve 



(2.10) 



(2.11) 



qvl -Ap^l ) \vi ) \ -i^tvo + 2p$^t;o,^ + p^^^v^ 
and find the expansion coefficients to be 

2/$xt^o,x , ^fx^x ^fgx^x , 2i<l>t i^t , 3/$^^ 

Ul= 2 + ^ T~2 + n + n + ' 

gvQ ^gvo ^g^vo 3gvo 3qvo gvo 

making use of the facts that uq = —2$^^ and p = q^. For j = 2 we solve 

-4/$^ gul W ^2 \ ^ / -huo - gufvo - 2guoUiVi - iuo,t - fuo,xx 
qvl -^p^l J \V2 ) \ -rvo - qvfuQ - 2qvoViUi + ivo^t - pvo,xx 

(2.12) 

and find the expansion coefficients U2,V2- However, we do not reproduce them here 
as they are so lengthy. 

At the resonance level j = 3, we have a linear relation between u^, and one of 
them is arbitrary. We see that the system 



2/$^ gut, \ u,\_ F, 



qvl -2p^i J J \G: 



(2.13) 



gives the compatibihty condition at j = 3 as 



qvoFs = guoGs. (2.14) 

Since the compatibihty condition has to be satisfied for any choice of the manifold 
coefficients of $t, ^tt occurring in fl2.14p must vanish identically. From the 

coefficient of ^tt we find that = q^. It follows that g = q or g = —q. The first is 
possible when g = gi and for the second we have g = ig2- Since / and g are both real 
or pure imaginary, we have f = fi for the former and / = 2/2 for the latter. 
We consider two cases: 

(i) The case g = ig2, f = if 2- From the coefficient of we find that 

2V0,x . g2 ,X J I.X n^XX r\ in ir\ 

^^0 g2 12 ^x 

which is a restriction on vq that contradicts its arbitrariness. We conclude that when 
/ and g are pure imaginary, the equation cannot pass the Painleve test. 

(ii) The case g = gi, f = fi- From now on, we are going to drop the subscript 1 and 
assume that / and g are real functions. Equation (12.141) is satisfied if and only if 

fA\f'-^'^Ji,V_f^_,i^^O. (2.16a) 



// fg \9 J f 9 

I. I r,! 9x r, 7 I 3 ftfx 1 fx9t , ftgx r,9t9x fxt , gxt ^ 

h2— + 2/i2 3/i2,x + -—PT ~ n-r- + 2^ H = 0. (2.16b) 

/ 9 2 2 f9 f9 9^ f 9 

Next we proceed to the resonance level j = 4. From the system 



2/$^ gu', ] u,\^ F, 
qvl 2p^l )\v, \G. 
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(2.17) 



we see that one of and should be arbitrary and it yields the compatibility 
condition 

voFi + uoGi = Q. (2.18) 

This relation contains the functions M3, ^3 of the resonance j = 3. From (12.131) we 
substitute U3 = {gulv^ — F3)/(2/$^) in (I2.18p . Similar to the previous resonance, 
from the coefficient of ^tt we get 



which integrates to give 



1^ + 2— = 0, (2.19) 
/ 9 



f{x,t)g\x,t) = K{t), (2.20) 



where K{t) is an arbitrary integration constant. Notice that this condition is 
already satisfied for the special case when / and g are only functions of t. We set 
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f{x,t) = K (t) / g"^ {x , t) and under this constraint see that fl2.16ap is satisfied whereas 
(]2.16bl) simphfies as h2^x = {dt/ g)x and determines the complex part of the potential 

/i2(x,t) = ^ + 7(t), (2.21) 
9 

where 7()f:) is an arbitrary function. Notice that /12 may depend on x if 5^ depends on 
X. The only remaining equation that has to be satisfied is (12.180 : 

9 \ 9 9^ ) " + 'k^^K \kJ ^Kg g ^ ^^^^^^ 

+ f-364 + 48^ - 6% - 10^^ + 

V ^ ^ f ^ ^ 

It is possible to integrate this expression and find real part of the potential as 
hi{x,t)= ^(^jgdx^ + j g^j r,dx)dx + ^ -^^]+6o [ gdx + 6i, (2.23) 



where Tq and ri denote 

272 {^K)t ^ {In K)u ... 1 fgt\ 

Here ^^{t) and 5i(t) are arbitrary functions. Note that integrals in fl2.23p will be 
evaluated without additional integration constants. 

Proposition 1 Equation (11.11) passes the Painleve test if and only if the conditions 
fl2:20|) . dm]) and (l2:23|) hold. 

We now consider a somewhat simpler case. We see from (12.191) that if one of / and g 
is only a function of t, then so is the other. Furthermore, we do not have a condition 
like (I2.20p for / and g, they can be independent. For / = /(t), g = g(t), ( I2.16ap 
holds identically and from ( ]2.16bl) we have /i2,x = 0, which means 

h2{x,t) = I3{t). (2.25) 

Similarly, we are left with a single equation from (l2.18p : 

2//1M. = 4/32 + 2/3 + 2/3^ - 4/3- - (7)' - 7- + 2(^)2 + ^ - ^. (2.26) 

/ 9 119 9 19 

This relation also follows directly from (I2.22p with the substitution 7(t) = /3(t) — |. 
We can write the right hand side as 

2//11,.. = 4/32 + 2/3 + 2/3^ - 4/3^ + (ln/)« + g{-)u - (In /)i(ln (?)i, (2.27) 

19 9 



where the subscript t denotes the derivative, and integrate to obtain hi{x,t): 

(2.28) 

with Hi{t) and H2{t) arbitrary. This conclusion thoroughly coincides with those of 
[IHl [23]. For / = 1 and ^ = -2 we have 

hi{x, t) = ^(2/3^ + /3)a;2 + Hi{t)x + H2it) (2.29) 

and it is compatible with the results of [2^ . 



2.2 Painleve test for eq. ( 11.21 ) 

We emphasize that there is no extra difficulty in applying Painleve test to other 
variants of (11. ip . One of the physically important possible extensions is (II. 2p . 
where /, g, h, k are as before complex-valued functions. For the special cases of the 
coefficients the equation is referred to as NLS equation with effective mass. If the 
singularity analysis is repeated for ( 11. 2p . it can be shown that passing the P-test is 
only possible when functions /, g are real and the following conditions hold: 

ki{x, t) =l-Ux, t) + fix, tf-^p^ (2.30) 
2 g{x,t) 

hix, t) =i (^(In + hiln + 2k2,x^ + -fit), (2.31) 

and 

H J \ 2^ K2,a; H J '^2,xt H K2K2,xx 

fhl ^ ^ffxQl _ ^Pgt ^ fxQxfxx _ ^ffxxQl ^ 'iflQxx ^"iffxQxQx 



(2.32) 



2^2 2g^ g^ 4g 2g^ 4g 2g^ 

9x9xx , '^ffxxgxx 2/ g^^ fgxfxxx , <^ffxgxxx 3/ gxgxxx , / 9x 



_|_ . JxUxx _|_ Jxxjjxx _ iJxx _|_ . ■^■^■^ _|_ . . ^^^^^ _ , ^^^^^^ _|_ 

9^ 9 9^ '^9 9 9^ 9 

Here, kxix, t) and k2ix, t) are real and imaginary parts of A;(x, t) and 7(t) is arbitrary. 
Two successive integrations of ( 12.32^ lead to the real part of the potential: 



hiix,t) = f ^ I Fsix,t)dxdx + 

J \/\.f\ J 



fx9: 



V0 r^''^*^"" 
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Here A{t) and B{t) are arbitrary functions and 



where e = 1 for / > and e = —1 for / < 0. 

We finally remark that at the time this work has been finalized it came to our 
attention that the recent paper [21] treated the variant of this equation in which / is 
real. 

2.3 Painleve test for eq. ( 11.31 ) 

We assume that the coefficients l{t),'m{t),n{t) are not identically zero. Furthermore, 
leading order analysis shows that 3m{t) + 2n{t) cannot vanish identically, since 
otherwise we would have l(t) = 0. We obtain integer resonances from the Painleve 
series expansion if (i) n(t) = —m{t) or (ii) n{t) = — |m(t). 

(i) We find the resonances as j = { — 1,0,1,3,4,5}. Compatibility conditions at 
the resonances hold if 



37 

/i(a;,t) = -^7i, /2(a:,t) = 0, g2{x,t) = (2.35) 
m 

31 31 
ki{x,t) = — gi^x , k2{x,t) = 5- 5'i + a;7i + 72 (2.36) 

h{x,t)= [ (^-!!^)dx+{xj, + ^2)--—,g! + -gi,.. + x/3,+/32 (2.37) 
J \ m / m m'^ m 

, . . rh I 31 ,„ 

where jijPiji = 1,2 are arbitrary functions of time. There is no special condition on 
gi{x,t),m{t) and l{t). 

(ii) Resonances occur at j = { — 1,0,2,3,4,4}. We obtain the same conditions of 
(i) with /3i(t) = in (K37\i . 

We stress that these results generalize those found in [22] in a consistent manner. 
For the case (i), considering the above results with 71 (t) = Pi{t) = and neglecting 
the indefinite integral in fl2.37p . one obtains exactly the same conditions as in the 
presence of time dependent coefficients. The same arguments lead to the conditions 
for (ii) in the time-dependent case, with the additional equation gi = Som, where 
So 7^ is an arbitrary constant. 
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3 Transformation to the Standard NLS Equation 
by Allowed Transformations 

In [20] the authors give a classification of the symmetry algebras of (11. ip into canonical 
forms using the allowed (or equivalence) transformations of the form 

^(x, t) = Q{x, t)i>{x, i), x = X{x, t), i= T{t), (3.1) 

which leaves form of (11. ip invariant, that is, the transformation does not add any term 
to the equation but may change the coefficient functions. Here Q{x,t) is a complex 
function obeying the constraint 

iQXt + f{x,t){XxxQ + 2Q,Xx) = 0. (3.2) 

If we set Q{x,t) = R{x,t)e^^^^'*\ R > 0,0 < 9 < 2n, the coefficient functions f,g and 
h map into 

mi) = f{x,t)^, 9{x,i) = 9{x,t)^^, (3.3) 



h{x, i) = l{h{x, t)-et + fix, t) - el 

+ 1 h2{x,t) + ^ + f{x,t){'^^ + e. 



(3.4) 



} 



through this transformation. Note that here f{x,t) is a complex- valued function. 
However, when writing h we have not separated it into real and imaginary parts, 
since from now on it is going to be taken as a real function. 

One of the main results of [20] is that any equation of the form (II. ip with a 
five-dimensional symmetry algebra is equivalent to the NLS equation with / = 1, 
g = e + ig2, h = 0, where e = ±1 and g2 = constant. The algebra is solvable and has 
a basis 

Po = di, Pi = d^, W = dc,, B = id^ + ^S:d^, D = idt + ^id,-^pdp, (3.5) 

which is isomorphic to the one- dimensional extended Galilei similitude algebra gs(l). 
Here ip is expressed in terms of the modulus and the phase of the wave function: 
ip{x,t) = p(x,t)e^'^(^'*). 

In the previous section we have been able to give the conditions for the equation to 
have the P-property. Now we ask, under these conditions, whether we can transform 
(II. ip to the standard NLS equation which is an integrable one from a group-theoretical 
point of view. This is equivalent to the condition that the equation under study has 
a five- dimensional symmetry algebra which is isomorphic to that of the standard 
NLS equation. Equivalence (or allowed) transformations will be used to produce the 
transformations mapping (II. ip to the NLS equation. 
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In Section 1, from the Painleve test it turned out that / and g have to be real 
functions. When / is a real function, (13. 2 p is equivalent to 



R\x,t) = ^, e. = -^^, (3.6) 

where -Ro(^) arbitrary. For the equation (11. ip to have a five-dimensional symmetry 
algebra, we only need to set 

/(x,t) = /(x,t)5 = l, (3.7a) 

~g{xrt) = 9{x,t)^k^ = e = ±l, (3.7b) 

/i(x,t) = 0. (3.7c) 

Notice that we have not included the constant ig2 in the righthand side of (I3.7bp since 
the functions occurring on the left are real. These equations are going to provide us 
the transformation functions i?, 9, X and T. 

We recall that we have succeeded in writing the potential h in terms of / and g 
for two (not distinct) cases. 

(i) The case when f{x,t)g'^{x,t) = K(t). 

In the following calculations we used the substitution / — )■ K/g"^. First solving 
()3.7ap we find that 



X, = ei]J^g, X{x,t) = ei]J^ j gdx + ^t), ei = ±1 (3.8) 



with an arbitrary ^. From the second equality ( I3.7bp we easily find 



• \ 1/2 

R{x,t) = I j • (3-9) 

Next we make use of the equations in (13.61) . The first equation means that the product 
R'^{x, t)Xx{x, t) should be a function of t, which is indeed the case. Integration of the 
second equation determines 9: 



9{x, t) = ^In ^ j ( j gdxf - ^ j dij gtdx)dx - ^^4^ j ddx + r/(t), 

(3.10) 
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where t] is an arbitrary function. The remaining equation to be solved is fl3.7cp . We 
are going to make use of fl2.2ip . fl2.23p . (13. 9p . (13. lOp in the formula (13. 4p . From the real 
and imaginary part of f]3.7cp we have, respectively, 



"iW + 1 + 1 = 0- (3-12) 

The second equation is easy to integrate and at once yields the time transformation 

r ^-'if'rd.t 

T{t) = Ti / ———dt + T2, Ti, T2 constants. (3.13) 
J K[t) 

Once /, g and h are given by a specific case of (11. ip , one can first check the P-property 
conditions and after that determine the function T, in principle. With this form of T 
the last term in (13. lip vanishes, but the first and the second ones do not. However, 
we can set them equal to zero by choosing ^ and i] as solutions of the equations 

W + ^iL--^^ = 0, (3.14a) 
2\/KT 2\/KT^I^ 

m-n-^ = ^- (3.14b) 

Having found T we can arrange fl3.14ap as 

of which integration twice yields 

m = -2ei JtiJ ^5oit)dt)dt + eor(t) + ei, (3.16) 
where ^1 are constants. From fl3.14bp follows directly 



r]it) = J 5r{t)dt- jf{j ^5,dtfdt + t,i, jf{j ^5,dt)dt-^T + e,, (3.17) 
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with 60 constant. So far we have all the information to write the general 
transformation formula. T and R are given by fl3.13p and (13. 9p . respectively. For 
the space variable transformation we have 

X{x,t) = j gdx-2e, j f{j ^5^{t)dt)dt + i^T{t)+^^, ei = ±1. (3.18) 

The phase of Q is obtained in the form 

^'^ " ~ 8^ ^] ^^^^^ j ^^j 



- jTiJ ^J^6odtfdt + ei^o Jt{J ^5odt)dt - |t + 9o. 

Now we make use of the formula we have found for T and give the transformation 
functions together with (I3.13P in its most explicit form. With T(t) = 2 J ■j{t)dt we 
have 



X{x,t) =tx^i— / gdx-2tx^x / ^( / e^KSodt)dt 

^dt + ^,, (3.20) 

i?(a;,t)=e-r(^-^j , (3.21) 

e{x,t)=^{2^+^){j gdxY-^ j g{j gtdx)dx- j ^-^{ J KSodt)'dt 

+ ^(| e'KSodt - ^i^) Jgdx + J 6,dt - ^ I '-^dt 

+ ei^oVTi J^iJ e^K6odt)dt + Oo- (3.22) 

(ii) The case when f{x,t) = f(t), g{x,t) = g(t). 

Although this case is already included in (i) we intend to present the results 
since a wide literature is devoted to the study of equations with this special form 
of the coefficients. In this case, eqs. (I3.14p are not imposed but rather satisfied 
identically. Transformations can be directly obtained from ( 13. 13p .( 131201) - (13. 22p via 
the substitution 

Kit) = fit)g\t), ^(t) = m-j^y Soit) = ^, 5,{t) = H,{t). (3.23) 
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Allowed transformations are found to be, with B{t) = 2 J f3{t)dt, 

T{t) =Ti J ^e-^^'dt + T2, (3.24) 

r 2 p f 



X{x,t) =eiVriye-^x-2eiVTi / —e''"{ / '-e''H^dt)dt 

+ ^oTi j ^e-2^rft + ei, (3.25) 
R{x,t) =e-^ (eTi^ , (3.26) 



/ 



9{x,t) =i-(2/3 + (ln^),)x^ + §e-^( / H,dt - ^i%^)a; 
4/ 9 f J 9 2 



eoTl f 9^23^^ 



^e-2^(/ LeBH.dtfdt- ^ , ^ 

f J 9 4: J f 

+ ei^oVT\ j J^'^^'i j ^e''Hidt)dt + J H^dt + 9^. (3.27) 

Remark 1 For f = fi + if 2, one can set = T / fi and this normalizes f through 
(13. 3p into f{x,t) = 1 + if2{x,t) 120^ . Thus, without losing any generality, we could 
equivalently investigate the equation 

i^t + (1 + if2)^.. + 9{^. t) 1^1 V + h{x, t)tp = 0. (3.28) 

From Painleve analysis of Section 2 would immediately follow /2 = and g = g{t) as 
a real-valued function and the potential components 



h,{x, t) = {P' + ^-P^ + ^{-)u)x^ + H,{t)x + H^it), h^ix, t) = Pit). (3.29) 
2 9 4: g 

Transformation formulae would considerably simplify in this setting. Yet, our results 
will remain useful when one begins with an equation with / 7^ 1. 



4 Integrability of a Gross-Pitaevskii Equation 

We apply our results obtained in previous sections to the equation 

iA + i^xx + g{x, t) + kx^i) = 0, (4.1) 

where g ^ C and G M. This equation governs the dynamics of a Bose-Einstein 
condensate and its integrability was studied in [T7j through the Painleve analysis. 
The equation was shown to pass the Painleve test for PDEs (WTC test) when the 
coefficient g satisfies 

gg-2g^ + Akg'' = (4.2) 
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and has the special form 



where A,B,go are arbitrary constants. Subsequently, by a proper transformation the 
equation is transformed to standard integrable nonlinear Schrodinger equation. 

We can recover the equation f l4.2p using our results of Section 1. Then, without 
considering results of the Painleve analysis, we will try to transform the equation into 
the standard NLS equation based on its group-theoretical properties. To this end, we 
shall require that the equation under study has a five-dimensional symmetry algebra. 
We then show that this invariance requirement will force g to have exactly the same 
form obtained by the Painleve approach. 

Identifying (14. ip with (II. ip . we see that for the Gross-Pitaevskii equation fi = 
1? /2 = 0, coefficient of the cubic term is the same complex function g{x,t) and 
hi{x,t) = kx"^, h2 = 0. Since / is a constant (which is equal to 1), g is necessarily 
a function of t only. We make use of (I2.26P and obtain the same equation (14. 2 p 
for g, under which equation (14. ip does have the P-property. If g is taken as a 
solution to (14. 2p . the GP equation transforms to standard NLS equation through 
the transformation formulae (I3.24p -( l3.27p . 

Existence of a transformation into standard NLS equation naturally motivates the 
question whether (14. ip can have a five- dimensional symmetry algebra for some form 
of coefficient function g{x, t), in particular the one given in (14. 3p . We try to transform 
(14. ip directly to the NLS equation via (13. 6 p and (13. 7p . From (I3.7ap we have 

X{x,t) = eiVfx + ^{t) (4.4) 
and (13. 6p implies R{x,t) = R(t) and 

e{x,t) = -llx'~'-^x + vit). (4.5) 
2VT 



(I3.7cl) has real and imaginary parts, both to be set to zero. From the imaginary part 
we have -|- O^^ = 0, which integrates to give 



R{t) = /^o[T(t)]^/^ Ro = constant. (4.6) 

From the real part we have d'^ + 6t — kx'^ = and this gives differential conditions on 
T(t),e(t) and r]{ty. 



.. \ 2 

T 3 / T 



8A; = 0, (4.7a) 



T 2 \ T 

^f-^T = 0, (4.7b) 
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i'^ + Afif = 0. (4.7c) 

The second equation yields 

e(t)=eoT(t)+ei (4.8) 
with ^0)Ci constants, while from the third it follows that 

with the constant ^o- These determine 6 in terms of T: 

9{x,t) = - '-fTnt)x- f T(t) + ^0 (4.9) 

There remains the solution of fl4.7ap . which is recognizable to be a Schwarzian 
differential equation which we rewrite as 

{T,t} + 8k = 0, (4.10) 

where {T, t} denotes the Schwarzian derivative of T with respect to the variable t. By 
the Schwarz theorem it is well-known that its solution is given by a ratio T(t) = 
where cxi and (T2 are linearly independent solutions of 

a(t) - 4A;a(t) = 0. (4.11) 



For a derivation of this result, the reader is referred to An alternative formulation 
based on a unimodular group invariance can be found in [22]. Furthermore, if the 
constants Ti, ...,T4 satisfy T2T3 — T1T4 7^ 0, Tiai{t) + T20"2(t) and T^ai^t) + Tia2{t) 
are also linearly independent. Hence 

is a general solution to f l4.10p . Linearly independent solutions of (14. lip have constant 
Wronskian, so, without loss of generality, we can arrange to have 

iy[Tiai(t)+T2a2(t),T3ai(t)+T4a2(t)] = {T^n- T2T^)W[ai{t),a2{t)] = -1, (4.13) 

which provides T as a solution with three arbitrary constants to the third-order 
equation fl4.10p . After further calculations, we are going to see that this choice of the 
Wronskian will enable us to obtain the function g{x,t) in its simplest form. 

Let A; > 0. From (14. lip we have (Ji{t) = cosh.2\/kt and o"2(t) = smh2\/kt and 
fl4.13p imposes T2T3 — T1T4 = We write 

^ Ticosh2v/fct + r2sinh2v^t 
Tscosh2y/kt + T4smh2y/kt 
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and 

f{t) = [Tscosh2\/kt + nsmh2\/kt]-^. (4.15) 

For < 0, we choose solutions of (14. lip as cri(t) = cos 2\/—kt and o'2{t) = 
sin 2\/—kt and write 

y^^^ ^ Ti cos 2v/^t + T2 sin 2y/^t ^^^^^^ 
T3 cos 2\/—kt + T4 sin 2^ —kt 

and 

T(t) = [T3Cos2v^t + T4sin2y^t]-2, (4.17) 

with the similar condition T2T3 — T1T4 = ^^^^^ . 

Now we have determined all the unknown functions of the equivalence 
transformation. Our aim was to determine g[x^ t) so that (14.11) has a five-dimensional 
symmetry algebra. From (14. 6 p and ( I3.7bp we find 



g{x,t) = git) = ^f'/'. (4.18) 
The following summarizes possible forms oi g{x,t): 



g^{t) = ^[Tscosh2Vkt + Tismh2y/kt]-\ k > 0, (4.19) 
R, 



■0 



e 



^2(t) = -2 [T3Cos2V-H + T4sin2V-H]-\ A; < 0. (4.20) 

Rq 

Now we would like to compare the results obtained by symmetry consideration with 
those given in [17]. In a search for the function g{x,t) through Painleve analysis, 
authors of [17] obtained the differential equation (14. 2 p and write (14. 3 p as its solution. 
Setting 

f = Lg^, L = constant (4.21) 

taking the relation (I4.18P between T and g into account, and plugging it in (I4.7ap . we 
obtain exactly the equation (14.20 for g. This implies that (14. 3 p ought to be contained 
in g{x,t) obtained by our approach. Indeed, if we choose 

e 2go e 2go 

in ( I4.19p . we see that it transforms to (14. 3p . We note that the solution ( 14. 3p does not 
include the case when k < 0. The possibility for g given by (I4.20p should be added 
as an integrable case. 

In fact, (I4.7ap is nothing but a special case of the third term which vanished 
identically in (13. lip because at the beginning of Section 3 we assumed that the 
equation have the P-property; i.e., (I2.23P holds. But now, since we have not taken the 
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P-property into consideration, it is obtained as an equation that has to be satisfied. 
We show that it is equivalent to the condition that the equation has the P-property. 
To illustrate our results in a tidy form we introduce the functions 

_ (Ti cosh2v^t + T2 smh2y/kt, for A; > 
[Ticos2v^t + r2sin2v^t, for < 0, 

J T3 cosh2v^t + T4 sinh 2A/^t, for /c > 
[T3Cos2v/^t + Ttsin2v/^t, for A; < 0. 

In conclusion, we have shown that, if g{x,t) of (14. ip is one of gjit), j = 1,2, using 
the transformation 

m = m = ^, 



x{x, t) = X(x, t) = ^x + + 6, 



R{x,t) 



Rq 



(4.24) 



we can convert (14. ip to the equation 

#i + fe + e|^|'Vi = 0. (4.25) 
Thus, if ip{x,i) is any solution of (14.251) . then for Qjit), j = 1,2 

t/j{x, t) = R{x, t)e^^(^'*) 4) {X{x, t),T{t)) (4.26) 

is a solution of (14. ip . 

As a by-product, we can transform the symmetry algebra (13.51) by the equivalence 
transformations to obtain a basis for the symmetry algebra of (14. ip as follows: 

= V% + {-^oV + VVx)d,-^pVVd,+ (^^-^Vx + ^iV^ + VV)x^^d^, 

X2 = Vd, + ^{-^o + Vx)d^, 

X3 = d^, (4.27) 

X4 = Ud, + ^{^^ + ^^^x)d^, 

X, = UVdt+(^-jU + ^V+{^ + UV)x^d,~^p{l + UV)dp 

+ f-^ + (-^1 + ^1/ - ^^). + 1(^ + ^ + UV)x^] 0.. 
' 4 ^ AV 4 4 V 4V V > \ ^ 
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5 Summary 



We took the approach of Painleve test combined with the symmetry properties as 
a prehminary attempt towards integrabihty of the most general variable coefficient 
NLS equation. The other tool which plays a central role in the analysis is the notion 
of equivalence transformations relating the VCNLS equation to its canonical form. 
Through the Painleve test we obtained restrictions on coefficients of the equation 
under study. This motivated us to construct suitable point transformations taking 
our original equation to standard cubic NLS equation. As confirmed by other similar 
studies existing in literature, equations passing Painleve test possess a Lax pair which 
is usually considered as an essential characteristic of complete integrabihty. We see 
that this is also the case for the coefficients singled out by the Painleve test. 

The results of this paper can be used to obtain exact (soliton) solutions of the 
VCNLS equation, when the coefficients have some specific form, from those of the 
integrable NLS equation by point transformations whose precise forms are presented 
throughout the paper. 
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